Diamagnetic Susceptibility and Current Distributions in Granular Superconductors at 

Percolation 
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We study a two-dimensional granular superconducting network at the percolation threshold under 
the influence of an external perpendicular magnetic field. By numerical simulations on the full 
nonlinear problem, we determine the scaling exponent for the magnetic susceptibility. Further, we 
report on the scaling properties of the current distribution. The scaling of the current is found to 
be independent of the value of the magnetic field. Our results are in contradiction with previous 
numerical results based on linearized equations. We find a value for the susceptibility exponent which 
does not agree with existing theoretical suggestions, but agrees perfectly with renormalization group 
calculations. 
PACS: 64.60.Ak,75.20.-g,74.50.+r,74.80.Bj 



Superconducting networks have been subject to inves- 
tigation in various topologies and approximations over 
the years. The granular superconducting network at the 
percolation threshold is the main concern of this text. In 
such a system under the influence of an external magnetic 
field, the scaling of the magnetic susceptibility obeys a 
power law. There exist two theoretical predictions for the 
value of the scaling exponent [1,2]. Previously numerical 
work has been done in order to decide on which theo- 
retical prediction is correct [3,4]. This work was approx- 
imate due to linearizations of the governing equations. 
We demonstrate in this paper that numerical work on the 
non-linear equations give different current distributions 
in the network. There exists further an estimate for the 
scaling exponent of the susceptibility based on renormal- 
ization group calculations [5]. Our simulations provide 
a value for this exponent that differs from the two theo- 
retical suggestions, but agrees with the renormalization 
group estimate. Attemps to compare with experimental 
result have been made, but they were not conclusive in 
favor of any theoretical prediction [6,7]. 

Bulk superconductors as well as granular supercon- 
ducting networks have been modeled by division into 
discrete nodes between which current can flow. Below 
a critical temperature T c the amplitude of the supercon- 
ducting order parameter is constant over the network, 
whereas the phase of the order parameter obtains, in 
principle, a different value at each node. The current 
flowing between two neighboring nodes is then a func- 
tion of the phase difference between the nodes. 

More generally, disordered or granular superconduc- 
tors are studied in the same way - as a regular network 
of nodes. The obvious physical difference in topology 
between the bulk superconductor and e.g. the random 
superconductor insulator-mixture constitutes no mathe- 
matical difference. Randomness and disorder are incor- 
porated by allowing each pair of neighboring nodes to be 
either connected or disconnected. The individual connec- 
tion between nodes, or the network contacts, are treated 
as Josephson junctions. Josephson junctions appear in 



many shapes. In this paper we use the so-called Resis- 
tively Shunted Junction(RSJ) model [8]. Some authors 
use the term RSJ when referring to the more general 
Resistively and Capacitively Shunted Junction(RCSJ) 
model. However, these models are fully equivalent in 
the physical situation considered here. 

In this paper we study the random superconductor- 
insulator in two dimensions at the percolation threshold 
[9] in a constant perpendicular external magnetic field. 
In equilibrium there will flow electric current in every 
closed superconducting loop in the network. This con- 
stitutes a critical phenomenon situation where the trans- 
port properties become highly nontrivial. Furthermore 
we emphasis the importance of the nonlinearity of the 
transport equations. In particular, we find the scaling 
exponents belonging to the moments of the current dis- 
tribution, and, more important, the scaling exponent of 
the diamagnetic susceptibility. 

The model is a two-dimensional square array of nodes. 
Each node is with probability p — 0.5 connected to each 
of its four nearest neighbors by bonds. All bonds are 
modeled as RSJ contacts. The phase of the supercon- 
ducting order parameter at node j is denoted <f>j. The 
constant perpendicular magnetic field B is represented 
by the vector potential A. We choose Landau gauge for 
convenience. Then the current ijk between two neigh- 
boring nodes j and k is 



ijk = ejk sin((/) J - - 4> k - Ajk) 



where 
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is the line integral from node j to node k of the vector 
potential. Here ejk is one or zero dependent on whether 
there is an RSJ contact between nodes j and k or not. 
The transport equations for the current in the network 
result from the Kirchhoff current law, for node j, 
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The sum is taken over the neighbor nodes. These equa- 
tions are equivalent to the ones obtained by minimizing 
the Hamiltonian of the system [10,4]. 

The diamagnctic susceptibility is given by [4,11] 
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where 



4>jk = 4>j - 4>k - A 
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Here the sum is to be taken over all nearest-neighbor 
pairs of nodes. S represents the area of the network. The 
susceptibility is obtained from (4) in the limit B — ► 0. In 
this limit the second term of (4) vanishes because ipjk is 
odd in B. 

Given a value for the magnetic field, the transport 
equations (3) can be solved for (fij , the phase of the or- 
der parameter at each node j. In turn the current in 
each bond is found from (1). The current distribution 
thus found is not multifractal. This is surprising as pre- 
vious work on the linearized case showed multifractality 
[12,13,9]. The concept of multifractals can be explained 
as follows. Having calculated the currents in the network 
bonds, the moments Mi of the current distribution are 
found, 



M i = S 4 = H e ^ [sin( 

<3-k> <j,k> 



cj> k -A jk )} 1 . (6) 



The sum is over all nearest neighbor pairs of nodes. Now, 
the question arises about how these moments scale with 
the linear size L of the networks considered. If each mo- 
ment obey a power law: 



M t ~ U 



(7) 



where the exponents r; do not have the form n = al + b, 
then the current distribution is multifractal. 

Hansen and Roux [4] have calculated the multifractal 
exponents resulting from the linearized version of (3), 
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It would seem that this approximation is valid for small 
B. We are interested in the behavior of the system in 
the limit L — > oo. This linearization process is in fact 
so that we first take the limit B — > and then the limit 
L — ► oo. That is not the same as taking L — > oo, then 
B — > 0. We report the results from the solution of the 
non-linear equations (3) and find a different set of expo- 
nents valid for all B. Likewise in the linear case, (8), a 
similar linearization of (4) and (6) shows that the dia- 
magnetic susceptibility is equal to the second moment of 



the current distribution. Whence the susceptibility scales 
as does the second current moment. 

This easy relationship between the current distribution 
and the susceptibility is a direct result of the lineariza- 
tion. There is no reason why this should hold also in 
the non-linearized case. Actually another set of coupled 
equations must be solved for (dipjk/dB) which appears 
in (4), before the susceptibility can be calculated. This 
set of equations is found by realizing that the left hand 
side of (3) is in fact an invariant with respect to the mag- 
netic field B. Thus the derivative of this expression with 
respect to B must be zero, for node j having k neighbors 
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By the solution of (3) the cosine factor in (9) is evaluated, 
and the equations become linear in the unknown partial 
derivatives. 

Numerically we have analyzed ensembles of seven dif- 
ferent system sizes. The networks in the ensembles are at 
the percolation threshold, meaning that each bond con- 
necting a neighboring pair of nodes is drawn by chance 
to be present with probability p = p c = 0.5. The equa- 
tions are solved using the conjugate gradient iterative 
scheme. The number of networks in the ensembles are as 
follows; L = 16 : 20000, L = 25 : 4000, L = 40 : 2000, 
L = 64 : 800, L = 100 : 130, L = 140 : 60, and 
L = 200 : 20. All the networks were analyzed for differ- 
ent values of the magnetic vector potential A. We show 
the calculated current moments in figure 1. Each of the 
curves correspond to one network size. It is intuitively 
clear that a larger network will carry more current than 
a smaller. On the other hand, that the exact scaling re- 
lation of Mi as a function of L should be a power law is 
not so clear. We also do not have much intuition for the 
value of the scaling exponent. The data in figure 1 can 
be rescaled with respect to their maximum value, shown 
in figure 2. 

The results show an excellent example of data collapse. 
This tells us that the scaling relation of the first current 
moment is the same for all values of the vector potential. 
The scaling exponent from this scaling is found from fig- 
ure 3. The upper curve corresponds to the scaling of the 
first current moment. The power law fit is very good, and 
the exponent is t\ — 2.01. The middle and lower curve 
in figure 3 correspond to the scaling of the second and 
fourth current moment respectively. Our work show that 
the scaling of all the current moments are independent 
of the value of the magnetic field. The scaling exponents 
found are summarized in table I. Here we also quote 
previous results based on the linearized equations [4]. 

These two sets of exponents differ drastically. The lin- 
earized equations give rise to a multifractal current dis- 
tribution. Our calculations on the non-linear equations 
show that those results must be wrong. The current dis- 



TABLE I. The table gives the scaling exponents of the mo- 
ments of the current distribution. The results are given for 
both previous linear results and present nonlinear results. 



Moment I 



n (linear) 



(non-linear) 



1 


2.26 


2 


3.03 


4 


4.75 


6 


6.68 


8 


8.61 



2.01 
2.00 
2.00 
2.01 
2.01 



tribution is not multifractal. It is not even fractal as the 
exponents within reasonable error limits are equal to the 
integer number two. 

The magnetic susceptibility % of the networks is shown 
in figure 4 as a function of the linear size L. A best power 
law fit gives the scaling exponent b/v = 0.90 ±0.03 where 
the error bar comes from fitting the data. In the same fig- 
ure we show with a dashed line the previously established 
result from the linearized numerical work, b/v = 1.03 
[3,4]. 

The previous result is in perfect agreement with the 
expression attributed to de Gennes [1], 

b/v = 2 - t/v = 1.025. (10) 

Rammal et al have given another suggestion [2], 

b/v = 2-t/v + (3/v= 1.129. (11) 

Here t = 1.300, v = 4/3, and j3 — 5/36 are exponents 
known from percolation theory [9]; conductivity critical 
exponent, correlation length critical exponent, and the 
critical exponent governing the vanishing of the density of 
the infinite cluster respectively. Wang and Lubensky has 
later calculated \ using a cumulant expansion method 
[5] . This renormalization group approach gives the value 
b/v = 0.91 ± 0.03, which is in perfect agreement with the 
present simulation. The experimentally quoted values 
are 0.99 and 1.09 [6,7], which are values comparable to 
all of the theoretical and numerical numbers. However, 
it seems not to be possible to draw conclusions on this 
basis. 

Our results are interesting for two major reasons. One 
being a demonstration of the dangers of linearizing a real 
nonlinear problem. This is equivalent to an interchange 
of limits as discussed above. Our results indicate that 
this interchange is not allowed, and the penalty is wrong 
answers. The other being that the value for the suscep- 
tibility critical exponent found is in contradiction with 
both existing theoretical predictions. The previous ver- 
dict was in favor of (10), but now the case must be re- 
opened. It might just well be the case that this exponent 
describes a fundamentally different property of the perco- 
lating system, other than the well known exponents v, t 
and p. Whence one should not expect it to be expressible 



FIG. 1. The first current moment as a function of the mag- 
netic vector potential Aq. We show the moments for different 
network sizes L as shown in the legend of figure 2. 



in turns of them. Generally it is known that renormaliza- 
tion group calculations can give reasonably precise esti- 
mates for critical exponents even though short cumulant 
expansions are used. At the same time standard series 
expansion methods of governing equations are known to 
fail in critical phenomenon situations. Therefore the re- 
sult of Wang and Lubensky seems more plausible, and 
the fact that numerical simulations gives the same result 
strongly indicates that this is indeed the correct value for 
the exponent. 

In summary, we have studied the granular supercon- 
ducting network at the percolation threshold under the 
influence of an external perpendicular magnetic field. We 
find the scaling exponent for the magnetic susceptibility 
to be b/v = 0.90 which differs both from existing theoret- 
ical suggestions and previous numerical work, but agrees 
with renormalization group calculations. Further, the 
moments of the current distribution is found to scale with 
exponent two, i.e., with the area of the network. This 
scaling is independent of the value of the magnetic field. 
This contradicts the previous studies of the linearized 
case which gives a multifractal current distribution. We 
conclude that this linearization is not legal and is equiv- 
alent to interchange of the limits B — > and L — ► oo. 
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FIG. 2. The same data as in figure 1, where each curve 
has been normalized with respect to its maximum value. We 
observe the good quality of the data-collapse. 



FIG. 3. Scaling of first, second, and fourth moments of the 
current distribution. The straight lines are best power law 
fits. The scaling exponents are found in table I. 



FIG. 4. The diamagnetic susceptibility \ i s shown as a 
function of the network size L. 
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